Universal and wide shear zones in granular bulk flow 
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We present experiments on slow granular flows in a modified (split-bottomed) Couette geometry 
in which wide and tunable shear zones are created away from the sidewalls. For increasing layer 
heights, the zones grow wider (apparently without bound) and evolve towards the inner cylinder 
according to a simple, particle-independent scaling law. After rescaling, the velocity profiles across 
the zones fall onto a universal master curve given by an error function. We study the shear zones 
also inside the material as function of both their local height and the total layer height. 

PACS numbers: 45.70.Mg, 45.70.-n, 83.50.Ax, 83.85.Cg 



Slowly sheared granular matter does not flow homo- 
geneously like a liquid. Instead, granulates form rigid, 
solid-like regions separated by narrow shear bands where 
the material yields and flows Q, S H S H @. Shear 
localization is ubiquitous in granular flow — think of ge- 
ological faults and soil fractures, avalanches, pipe flows 
and silo discharges 0,0,100000313. 

Despite their crucial importance, granular shear flows 
are still poorly understood, in part because shear local- 
ization itself remains enigmatic [II, 0, IE 9 ■ O n the one 
hand, shear bands have a typical thickness of five to ten 
grain diameters and such steep gradients are difficult to 
cap ture b y continuum theories |E IE IE IE IE IE 
EE IH III 13 El On the other hand, the experimental 
handles for probing shear localization are limited. For 
example, studies in Couette cells always show the for- 
mation of a narrow shear band near the inner cylinder, 
irrespective of dimensionality, driving rate or details of 
the geometry jj, IlU. [lJi [lj], Ha] : Shear banding is very 
robust. 

In this Letter, we introduce a general experimental pro- 
tocol that can yield very wide shear zones away from 
the side-walls. We modify a Couette cell by splitting its 
bottom at radius R s . The resulting concentric rings are 
attached to the stationary inner and rotating outer cylin- 
der, respectively, and the cell is filled with grains up to 
height H (Fig. When driving the system, a shear 
zone is found to propagate from the slip position R s to- 
wards the surface where we measure the average grain 
velocities. Note that our strategy differs from previous 
works, which were carried out for large filling heights and 
smooth bottoms so as to minimize the effect of the bot- 
tom boundary [E 0, EJ- Here we take advantage of 
gravity and " drive the system from the bottom" . 

Our main finding are: (i) For large H, a regime of 
wall-localized shear band near the inner cylinder is recov- 
ered [EEiEl, but for intermediate H, we observe shear 
zones of tunable width away from the boundaries. This 
paper focuses on describing these "bulk" shear zones, (ii) 
The angular velocity profiles ui(r) of bulk shear zones fall 
onto a universal master curve which is best fitted by an 
error function. These profiles are therefore fully charac- 
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FIG. 1: (a) Schematic side-view of our split-bottomed Cou- 
ette cell, showing the stationary bottom disk and inner cylin- 
der (both grey), the rotating bottom ring (striped) and outer 
cylinder. A thin layer of felt between the rings assures smooth 
rotation. The inner and outer cylinder radii are fixed at 65 
mm and 105 mm respectively, while R s can be varied, (b) 
Setup in disk geometry, (c) Top view of the filled set up in 
the Couette geometry, where the rectangle indicates the area 
recorded by the video camera. 



terized by two parameters only: their center position R c 
and width W. A concise presentation of these results 
has appeared in ^E|- (Hi) R c and W depend on H, R s 
and particle properties in specific manners. The center of 
the shear zone R c evolves to the inner cylinder with in- 
creasing H in a particle independent manner. The shear 
zone width W grows continuously with H and depends 
on particle size and shape, but not on the slip radius R s . 
(iv) For a given height inside the material, the width and 
position of the shear zones depend on the height of the 
free surface H . 

Setup - A sketch of our split-bottomed Couette cell is 
shown in Fig. When the inner cylinder of the "Cou- 
ette" geometry is removed we obtain the " disk" geometry 
(Fig. njj) . Different sets of bottom rings allow us to vary 
R s from 45 mm to 95 mm. Grains, similar to those used 
in the bulk, are glued to the side walls and bottom rings 
to obtain rough boundaries. We studied spherical glass 
beads of size distributions 0.25 — 0.42 mm (I), 0.56 — 0.8 
mm (II), 1 — 1.2 mm (III) and 2 — 2.4 mm (IV), and irreg- 
ularly shaped plastic flakes (1.0—1.6 mm) (V), aluminum 
oxide beads (1.5 — 2 mm) (VI) and coarse sand (1.2 — 2.4 
mm) (VII). After filling the cell, an adjustable blade flat- 
tens the surface at the desired height. The outer cylinder 



2 




65 70 75 80 r ( mm ) 65 70 75 r (mm) 

FIG. 2: Main features of the normalized angular surface ve- 
locity in the Couette geometry for 0.3 mm glass beads (mix- 
ture I) and R s = 85 mm. (a) u)(r) for a range of equidistant 
heights h — 3,8, ... ,53 mm (right to left), (b) Contour plots 
of io(r), where the symbols correspond to, from left to right, 
u> — 0.1, 0.25, 0.5, 0.75 and 0.9. The curve indicates the strain 
rate maximum and shows the rapid qualitative change of the 
profiles when the inner cylinder is approached. 



and its co-moving ring are then rotated. A Pulnix TM- 
6710 8-bit CCD camera records 2000-frame movies of the 
resulting flow at the top surface at a rate of 120 frame/s 
with pixel resolution 100 /im. 

The flow rapidly (~ 1 s) reaches a stationary state 
where it is purely in the azimuthal direction, so that the 
surface velocities is a function of the radial coordinate 

only HE HE HiEi. We checked that these veloci " 
ties are proportional to the driving rate fl 0, 0, 0, 

for 0.16 < Q < 1.5 rad/s, and subsequently fix Q at 0.16 
rad/s. We thus focus on the velocity profile uj(r), the 
dimensionless ratio of the average angular velocity and 
O. We measure u>(r) with high radial resolution by par- 
ticle image velocimetry, i.e., by determining the averaged 
angular correlation function as function of r of two tem- 
porally separated frames. Unless noted otherwise, the 
time separation between frames is around 0.3s. 

Basic phenomenology - Figure [3 illustrates the main 
features of these velocity profiles. For shallow layers, a 
narrow shear zone develops above the split at R s . When 
H is increased, the shear zone broadens continuously and 
without any apparent bound. The broadest observed 
zones exceed 50 grain diameters in width. Additionally, 
with increasing H, the shear zone shifts away from R s 
towards the center of the shear cell [2jj. Indeed, for 
sufficiently large H, the shear zone reaches the inner 
wall, where it approaches the asymptotic reg ime of wall- 
localized shear bands reported earlier [ol lid Il7j . Before 
this wall localization occurs, however, there is a substan- 
tial range of layer heights where wide and symmetric bulk 
shear zones can be observed. 

Universal velocity profiles - Figure |31 illustrates our 
main result: After proper rescaling, all bulk profiles col- 
lapse on a universal curve which is extremely well fitted 
by an error function: 

u{r) = 1/2 + 1/2 erf {{r-R c )/W) . (1) 

A residue analysis comparing the fit to Eq. (J| to an 
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FIG. 3: (a) Collapse of all 35 bulk profiles obtained for H = 
4, 5, 6, . . . 38 mm (mixture I, R s = 85 mm, as in Fig. |2Jl when 
plotted as function of the dimensionless coordinate A := (r — 
R c )/W. The grey curve is an error function, (b) Differences 
£ between bulk profiles and fits to Eq. (c) Differences 

£ between bulk profiles and fits to 1/2 + l/2tanh(ZA), where 
I « 1.20274 compensates for the difference in scale between 
error function and hyperbolic tangent. The grey curve shows 
(erf(A) — tanh(ZA))/2. (d) Left tail of the velocity profile in 
log-lin scale, compared to the best fit to either a hyperbolic 
tangent and error function (grey curves). Dots, crosses and 
pluses are obtained for time-lags 0.3, 3, and 30 s respectively. 



alternative fit to an hyperbolic tangent shows that the fit 
to the error function is always better (Fig. 3b-c). By re- 
peating this procedure for the other particle mixtures, we 
establish the general superiority of Eq. Q : particle shape 
does not influence the functional form of the velocity pro- 
files. The robust form of ui{r) contrasts with the particle 
dependence found for wall localized shear bands |5|]. For 
these, the vicinity of the wall causes layering, in particu- 
lar for monodisperse mixtures. Apparently such layering 
effects play no role for our bulk shear zones. Accurate 
measurement of the tail of the velocity profile (Fig. 01) 
further validate Eq. and rule out an exponential tail 
of the velocity profile here. The strain rate is therefore 
Gaussian, and the shear zones are completely determined 
by their centers R c and widths W. 

What limits the universal regime? Apart from wall- 
localization (see Fig. [2), we find that in the disk geome- 
try ui(r) starts to deviate from Eq. JI} when H exceeds 
~ R s /2. The symmetry of the velocity profile, easily de- 
tectable by a simple x 2 test, is then weakly broken [2l| . 
In the following, we focus on the functional dependencies 
of R c and W on the parameters R s , H and particle type 
for the universal profiles given by Eq. ^ 

Shear zone position- Remarkably, the shear zone cen- 
ter evolution with height, R C (H), turns out to be inde- 
pendent of the grain properties (Fig.Hk)- Therefore, the 
only relevant length-scales for R c are the geometric scales 
H and R s . The dimensionless displacement of the shear 
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FIG. 4: Shear zone positions i? c versus height, where H 
is restricted to the universal regime detected by a \ 2 test, 
(a) Comparison between R C (H) for glass beads I-IV (closed 
symbols), non spherical particles V-VII (open symbols) and 
Eq. @ (curve). For all runs, R a is 95 mm. (b) R C {H) for 
mixture III and R 3 — 95, 85 and 75 mm in the Couette ge- 
ometry (open symbols) and R s = 95, 65 and 45 mm in the 
disk geometry (closed symbols), compared to Eq. 2 (curves). 



zone, (R s — R c )/R s , should thus be a function of the 
dimensionless height (H/R s ) only. The simple relation 



(R s - R c )/R s = {H/R s f' 2 



(2) 



fits the data well (Fig.Ek) H^. To check Eq. 101 we have 
varied R s over a substantial range. Only the presence 
of the inner cylinder limits the range of R s . We find no 
differences between bulk velocity profiles measured with 
or without the inner cylinder — Bulk shear zones are 
insensitive to the presence of the side walls. So, we sub- 
sequently switched to the disk geometry (Fig. and 
obtained an excellent agreement between R c and Eq. (J2J) 
over the range 45 < R s < 95 mm ( Fig. ^J)) . 

Shear zone width - The width of the shear zones de- 
pends on the particle size and type (Fig. |3^,-c), but not 
on R s (Fig. |SJi). First of all, W grows with H and in- 
creases for larger particles (Fig. |SJi-b) . The data shown 
in Fig. Et-b can be made to collapse when plotted as 
W/d vs H/d (not shown), where d denotes the grain size. 
Grain shape and type also influences W(H): irregular 
particles display narrower zones than spherical ones of 
similar diameter (Fig. |5j;). Finally, for the universal ve- 
locity profiles, W is independent of R s (Fig. 01) ■ We 
therefore conclude that the relevant length-scale for W 
is given by the grain properties. 

The evolution of the velocity profiles from a step func- 
tion at the bottom to an error function at the surface, is 
reminiscent of a diffusive process along the vertical axis. 
However, W grows faster than y/~H as diffusion would 
suggest, but slower than H. Intriguingly, we obtain the 
best fit for W ex H 2 / 3 over the limited range where we 
have reliable data. We cannot, however, rule out other 
functional dependencies such as a crossover from square- 
root to linear behavior. 

Below the surface — So far we have only discussed ob- 
servations of the surface flow. To get some insight into 
the 3D bulk flow structure, we put patterns of lines of 
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FIG. 5: Width of the shear zones versus H. (a) Linear 
and (b) log-log plots of W for spherical glass beads of in- 
creasing sizes for R s = 95 mm (Mixture I-IV). The width 
grows with particle size. The straight lines in (b) have slopes 
1/2, 2/3 and 1, showing that W grows faster than \f~H but 
slower than H . (c) Irregular particle shapes diminish W sub- 
stantially, as can be seen by comparing the widths for the 
plastic flakes (V/pluses), aluminum oxide beads (Vl/open cir- 
cles) and coarse sand (VH/crosses) to those of the glass beads 
shown in panel (a), (d) W for Mixture III and i? s =95, 85 and 
75 mm (Couette geometry; open symbols) and R s ~95, 65 and 
45 mm (disk geometry; closed symbols). The strong upward 
deviations observed in the disk geometry coincide with the 
symmetry breaking of the velocity profiles (see text). Apart 
from this, the setup geometry does not affect W. 



colored tracer particles at given Hi, (Fig. |SJi) . More ma- 
terial is carefully added so that the line-pattern is buried 
under a given amount of grains (H > Hi,). We then ro- 
tate the system for a short period (sim8 s), and recover 
the deformed line-pattern by carefully removing the up- 
per layers of grains (Fig. Os) . Comparing the snapshots 
of the deformed pattern to the initial one allows for the 
determination of the velocity profiles in the 3D bulk of 
the material (Fig. EJ:). We have checked that transient 
effects are limited, that the measurements reproduce well 
and that there is no significant motion in the vertical di- 
rection. 

The position of the shear zones in the 3D bulk are 
presented in Fig. Clearly, the evolution of R c with 
Hf, inside the material depends on the total amount of 
matter, as given by H: the more material is added, the 
more the shear zone shifts towards the center. This ob- 
servation is confirmed by recent theory |23l | , numerics [iij 
and MRI measurements |25j . The widths of the 3D shear 
zones are more difficult to measure accurately, but a clear 
trend can be identified: shear zones become wider when 
more matter is added on top (Fig. Eli) . 

Outlook — In this work we have presented a simple 
experimental protocol in which wide and tunable shear 
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FIG. 6: Measurements of u)(r) inside the material for the disk 
geometry and R s — 95 mm. (a) Initial pattern at Ht — 31 
mm. (b) The same pattern after adding material up to H — 
49 mm, rotating and removing the upper layers, (c) Velocity 
profiles for H = 49 mm, and Hi, = 7, 13, 19, . . . , 49 mm (from 
right to left), (d) R c {H b ) for H = 25 mm (pluses), H = 37 
mm (crosses) and H = 49 mm (circles) compared to R C (H) 
given by Eq. JS) (curve), (e) Corresponding W(Hb). 



zones can be generated in a variety of granular materials. 
Perhaps the biggest surprise is the robust and remarkably 
simple form of the velocity profiles — for granular sys- 
tems, universality is rare. Our measurements of the tail 
of the profile indicate that even far away from the shear 
zone, the grains are not entirely at rest. These findings 
indicate that these features are amenable to a continuum 
description of granular matter 0, 0, 0, , in particular 
for very wide shear zones. 

It is also noteworthy that the functions for the shear 
zone location R C (H) and width W(H) depend on entirely 
distinct sets of parameters: the relevant length-scales for 
R c and W appear to be well separated. The particle- 
dependence of the width provides a characteristic length- 
scale which may bridge microscopic and coarse-grained 
descriptions. The development of theories of granular 
flows can further be guided by the universal relation 
for the shear zone position, and should incorporate the 
strong influence of the boundary: Avoiding the proxim- 
ity of the side wall can turn the shear zones from narrow 
to wide, and from particle dependent to universal. 

The broad shear zones occuring in our geometry allow 
for further experiments that are more difficult to real- 
ize in narrowly localized shear bands. Important issues 
can thus be tackled, such as the velocity fluctuations and 
particle diffusion for various strain rates and locations 
within the shear zone. One can also probe whether "local 
clusters" of particles, possibly similar to those found in 
more rapid flows, would occur for these denser flows |26| . 
Finally, over the range studied, the velocity profiles are 
rate independent, but what happens for much larger and 
smaller rotation rates is an open question. The simple ex- 
perimental protocol that we provide for creating generic 



(i.e., away from sidewalls) shear zones can be the starting 
point for many crucial experiments, thus addressing the 
basic question: "How does sand flow?" . 
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